Collinear Four- Wave Mixing of Two-Component Matter Waves 
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We demonstrate atomic four- wave mixing of two- component matter waves in a collinear geome- 
try. Starting from a single-species Bose-Einstein condensate, seed and pump modes are prepared 
through microwave state transfer and state-selective Kapitza-Dirac diffraction. Four- wave mixing 
then populates the initially empty output modes. Simulations based on a coupled-mode expansion 
of the Gross-Pitaevskii equation are in very good agreement with the experimental data. We show 
that four- wave mixing can play an important role in studies of bosonic mixtures in optical lattices. 
Moreover our system should be of interest in the context of quantum atom optics. 

PACS numbers: 03.75.Mn, 03.75.Gg, 67.85.Hj 



Four-wave mixing is a fundamental, well-studied con- 
cept in nonlinear optics and spectroscopy pQ . Its matter- 
wave analogue, based on binary collisions in ultracold 
atomic gases, was first demonstrated experimentally a 
decade ago [2j|3], establishing the field of nonlinear atom 
optics [4 j. In four- wave mixing (FWM), two waves form a 
grating from which a third wave diffracts, thus generating 
a fourth wave. This process has been used for coherent 
matter- wave amplification [5 , and for the generation of 
correlated atom pairs [4- 7 . Energy and momentum con- 
servation require the magnitudes of all atomic momenta 
in the center-of-mass frame to be equal which, for atoms 
in a single internal state, necessitates a two-dimensional 
geometry [3]. By modifying the dispersion relation with 
an optical lattice, nondegenerate FWM of a single species 
becomes possible also in one dimension [8]. 

Despite considerable theoretical work on atomic FWM 
with more than one internal state [9HTT], experiments 
have only very recently started to explore possible mech- 
anisms for such FWM [7J [12]. The additional internal 
degree of freedom allows for degenerate FWM to occur 
in one dimension, with pairs of waves in different inter- 
nal states sharing the same momentum mode, opening up 
possibilities to generate nonclassical matter- wave states, 
e.g. with macroscopic spin entanglement [10]. In this Let- 
ter, we demonstrate free-space collinear atomic FWM in- 
volving two internal states with distinct, macroscopically 
populated momentum modes. 

Apart from the relevance for quantum atom optics, an- 
other important context arises in experimental studies of 
bosonic mixtures in optical lattices [13-15 . These sys- 
tems are of high interest not only in connection with 
applications to quantum magnetism [16 , but also for 
studies of decoherence mechanisms [17], and for lattice 
thermometry [14] [15] . Most experiments with ultracold 
atoms in optical lattices to date rely on time-of-flight 
information. In particular, a sudden release from the 
lattice projects the band populations onto plane-wave 
states [18]. We find that for a homonuclear mixture of 
interacting superfluids, FWM processes can alter the ex- 
pected momentum-space distributions, masking or even 



mimicking in situ interaction effects. 

In order to induce collinear two-component FWM, we 
apply a state-selective optical lattice pulse to a Bose- 
Einstein condensate containing atoms in two internal 
states | a) and |6), as illustrated in Fig. [I (A). The pulse 
induces Kapitza-Dirac (KD) diffraction 19] [20] produc- 
ing recoiling \a) atoms in both positive and negative mo- 
mentum modes |±2) = |±2ftfcc) where kL = 2tt/\l, 
while the \b) atoms remain unaffected. Subsequently, as 
illustrated in Fig.[l](B), the \b) atoms Bragg diffract from 
the density modulation formed by the interference of the 
recoiling \a) atoms \a, 2) = \a) (g) |2) with those at rest, 
|a,0). Because of momentum exchange collisions, the 
recoiling \a) atoms are coherently transferred back into 
|0), as recoiling atoms |b, 2) are produced. This process 
is formally not distinguishable from coherent (pseudo-) 
spin exchange. Our system might thus pose an inter- 
esting alternative to spinor condensates for the creation 
of nonclassical states [lOj [12] . We note that due to the 
symmetry of the KD pulse, another, independent "copy" 
of the FWM process occurs on the negative momentum 
side. For quantum atom optics purposes, this can easily 
be avoided by using a state-selective Bragg pulse instead, 
which also allows for extended control of the initial mode 
populations. In the present work, however, we are con- 
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FIG. 1. Experimental scheme. (A) State-selective Kapitza- 
Dirac diffraction of a two-component Bose-Einstein conden- 
sate. (B) Four-wave mixing (solid arrows) with pump modes 
\b, 0), |a, 2) and seed mode |a, 0) transfers \b) atoms to the 
output mode \b, 2). Because of the symmetry of the prob- 
lem, the process also occurs for the modes |a, —2) and |6, —2) 
(dashed arrows). 
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tent with applying a KD pulse, mainly out of technical 
convenience. 

Our experimental setup has been described in detail 
in Ref. [2T] . In a crossed-beam optical dipole trap at 
1064 nm wavelength, we produce nearly pure 87 Rb Bose- 
Einstein condensates in the \a) = \F = l,rap = —1) hy- 
perfine state typically containing about 1.6 x 10 5 atoms. 
The trap is approximately isotropic with a mean trap fre- 
quency around 50 Hz and an alignment-dependent verti- 
cal frequency uj z /2tt between 40 and 50 Hz. Immediately 
after a variable fraction of the condensate is transferred 
into the state \b) = |2, — 2) via a microwave Landau- 
Zener sweep [22 , a state-selective lattice beam [23] at 
Xl = 785.1 nm is pulsed on along the vertical z direction 
(1/e 2 radius 230 /im) for a time tkd- The polarization 
(<r + ) is chosen such that only the \a) atoms feel the op- 
tical lattice potential formed through retro-reflection of 
the beam. A magnetic field (~0.4 G) along the beam 
axis defines the quantization axis. After release from 
the trap and a few milliseconds of free evolution, dur- 
ing which the FWM occurs, a magnetic field gradient 
(Stern-Gerlach pulse) spatially separates the two hyper- 
fine states along the horizontal x axis for detection. The 
atoms are imaged after a total time of flight of 15 ms 
via near-resonant absorption imaging by a 100 /is long 
pulse of F = 2 — >• F' = 3 imaging light, combined with 
F = 1 — >• F' = 2 repumping light, which ensures equal 
detection efficiencies for both hyperfine states. 

In Fig. [2] typical absorption images are shown for three 
different fractions of \a) atoms f a = N a /N. The KD 
pulse duration (25 fis) and the lattice depth V a for atoms 
of type | a) (6 Er, where Er = h 2 k\/2m is the recoil en- 
ergy), are chosen such that half of the \a) population is 
diffracted into |±2), while higher orders remain largely 
unpopulated. By analyzing single-component diffraction 
patterns [20J, we have determined the lattice depths for 
each component, confirming that atoms of type \b) expe- 
rience < 5% of the lattice depth seen by the \a) atoms. 
On their own, the \b) atoms therefore are not affected 
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FIG. 2. Typical absorption images taken after application of 
the state-selective Kapitza-Dirac pulse (jkd — 25 /is, V a = 
6 Er), 15 ms time of flight and Stern-Gerlach separation 
along x, for the case of (A) only \a) atoms present (f a = 1), 
(B) only | b) atoms (f a = 0), and (C) equal populations of 
both components (f a — 0.5). 



by the lattice pulse, as shown in Fig. [2] (B). However, 
when both components are present, a significant frac- 
tion of | b) atoms is transferred into the |±2) momentum 
modes [Fig. [2] (C)]. 

We have measured the amount of diffracted atoms in 
each state as a function of f a . As shown in Fig. [3] (A), 
the fraction of diffracted \b) atoms + -^6,-2)/^, 

monotonically increases from zero towards a maximum 
as f a is increased, consistent with the picture that the 
grating formed by interference of the |a,0) and |a, ±2) 
modes, from which the atoms in |6, 0) diffract, gets deeper 
as the number of \a) atoms grows. The relative number 
of diffracted \a) atoms has a pronounced minimum near 
f a = 0.5, which can be interpreted as a "backaction" of 
the I b) atoms onto the \a) grating. 

To obtain a more quantitative understanding, we the- 
oretically model our system starting from the coupled 
Gross-Pitaevskii equations (GPE) for the order parame- 
ters 3>Qi(r, t) of the two components a G {a, b} 

\ m Pe{a,b} J 

where g a p = 47rh 2 a a ^/m^ m is the atomic mass, and the 
intra and interspecies s-wave scattering lengths a aa , a^, 
and a a b in units of ao are 100.4, 99.0, and 99.0, respec- 
tively [24]. The trapping and lattice potentials are given 
by V2 ot = T4 r apM) + V a (t) sm\k L z). Similar to the 
slowly varying envelope approximation (SVEA) [3j[TT], 
we approximate the solution of Eq. as an expansion 
in terms of momentum modes, or wave packets, moving 
along z with multiples of the recoil velocity vr = Kk^jm 

00 

Mr,t)= E c na (t)e ink ^^ (r-znv R t,t). (2) 

n=— oo 

We further assume that the wave packets $o are of 
Thomas-Fermi form and that they expand hydrodynam- 
ically after release from the trap [25] , which leads to a sig- 
nificant simplification compared to a full SVEA simula- 
tion. On the time scales of interest, phase-separation [26] 
can be neglected, and we have $> a ex <l>o for both com- 
ponents just after the microwave transfer [26]. Since the 
momentum spread of <3>o is much less than te^, the modes 
in the expansion are quasi-orthogonal. After inserting 
the ansatz ([2| into Eq. ([T]), we arrive at a system of cou- 
pled equations for the amplitudes a n (t) = c na (t) 

ihd t a n = E R n 2 a n + V a (t) [ \ a n - \ (a n+2 + a n _ 2 )] 

+ ^ (gaa ^rn a rn / + dab b^b m f)a n f h nrnm > n >(t) 
mm' n' 

and similarly for the other component b n (t) = c no (t). 
Here, h nrnrn > n >(t) ex S(n + m — m! — n') denotes overlap 
integrals that include the effective temporal decay of the 
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nonlinear interaction, as the different wave packets sepa- 
rate, and as the density decreases during the expansion. 
The terms responsible for FWM are of the form b^b n a m 
(and a^a n 6 m for the \b) component) with m ^ n. Af- 
ter adjacent modes {\m — n\ = 2), for which the overlap 
decays the slowest, have completely separated, the pop- 
ulations remain frozen, since only equal-momentum self 
and cross-phase modulation terms of the forms \a n \ 2 a n 
and \b n \ 2 a n survive. With a typical Thomas-Fermi ra- 
dius R z ~ 10 /im, we obtain a typical separation time 
t sep « 2R z /2vr of 1.7 ms. 

The full set of observed populations |a n | 2 and \b n \ 2 
after FWM is plotted in Fig.[3](B), along with predictions 
of our model obtained with parameters according to the 
experimental ones, leaving only the total atom number 
N as a fit parameter. The overall agreement between 
data and theory is remarkable. The maximum FWM 
yield occurs near f a = 2/3 where the initial populations 
of the pump and seed modes are equal, maximizing the 
FWM term a^a n 6 m at t = [3 . The data also clearly 
show the correlated growth of |&±2| 2 and |an| 2 , along with 
a corresponding depletion of the pump modes |6o| 2 and 
|a±2| 2 , as detailed in the inset. 

We note that since the FWM yield is proportional to 
the interspecies scattering length as well as to the overlap 
jdr |^> a | 2 |$5| 2 of the two components, it can serve as a 
sensitive probe for both quantities. As a practical exam- 
ple, we use two-component FWM as a clear "single-shot" 
diagnostic for the optimization of component overlap. By 
carefully canceling magnetic field gradients, we are able 
to sustain overlap, i.e. FWM yield, for up to 2 s after the 
microwave transfer. 

To further confirm the coherence of the observed two- 
component FWM as implied by our model, we directly 
map out the time evolution of the output mode popu- 
lation \b±2\ 2 by interrupting the FWM process after a 
variable time through the selective removal of \a) atoms 
with a 50 /is long "blast" pulse of repumping light. As 
shown in Fig. [3] (C), the atom number in the |b, ±2) 
modes smoothly grows from zero to a maximum value 
reached around the expected separation time. The non- 
linear, initially quadratic growth is indicative of a coher- 
ent process [3j [5] (other signs would be an overshoot and 
oscillations, which however would require higher densities 
or longer overlap). To exclude the possibility that the 
observed growth is merely an artifact caused by density- 
dependent losses of | b) atoms accompanying the blast 
(due to collisions with \a) atoms), we repeat the exper- 
iment with the polarization of the lattice beam chosen 
such that both components experience the same lattice 
depth of about 6 Er. In this case, we expect the |6, ±2) 
modes to be populated immediately after the KD pulse, 
as indeed can be seen in Fig.|3](C). Further, no FWM 
is expected to occur for Vb = V a , as the internal and ex- 
ternal state dynamics are decoupled. By comparing the 
observed time evolution for this reference case with the 
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FIG. 3. Mode populations after four- wave mixing. (A) Frac- 
tion of atoms with momenta ±2hkL in state \a) (green circles) 
and | b) (orange squares). (B) Populations |c nc *| 2 = N na /N of 
the modes |a,n), as indicated in the plaquette. The dashed 
lines (1), (2), and (3) indicate the initial conditions before 
FWM, |6 (0)| 2 , |a (0)| 2 = |a± 2 (0)| 2 , and |6± 2 (0)| 2 (where ±2 
indicates the combined populations) . The arrows indicate the 
temporal evolution of the populations. The solid lines repre- 
sent the predictions of the model (V a = 5.6 Er, tkd = 25 /xs, 
uj z = 2n x 51 Hz, N = 1.4 x 10 5 ). The inset shows the trans- 
ferred | b) population vs the transferred \a) population, where 
the dashed line represents a slope of unity. (C) Growth of 
the population in \b, ±2) following the Kapitza-Dirac pulse 
for Vb = (orange squares) and for Vb = V a (blue triangles, 
xl/2). The FWM was interrupted after a variable time r by 
blasting away the \a) atoms. Each data point is averaged over 
2-6 runs (here uo z — 2n x 41 Hz, f a = 0.5). The dashed lines 
are the predictions of the uncorrected model (including higher 
order FWM terms), whereas the solid lines take into account 
the loss of atoms during the blasting process (see text). The 
blast-loss model was calibrated by fitting to the Vb = V a data. 



expected one, we can calibrate our model for the blast- 
induced losses, which assumes a relative loss of \b) atoms 
proportional to the density of \a) atoms in the overlap re- 
gion. With this correction, the theoretical time evolution 
for Vb = matches the experimental data very well. 

So far, we have discussed controlled FWM after ap- 
plication of a short optical pulse to induce diffraction. 
Now, we turn to the question whether FWM is also rele- 
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FIG. 4. Four- wave mixing effects for an adiabatically ramped- 
up optical lattice. (A) An | a) -selective lattice is ramped up to 
a depth of 6.0 Er within 100 ms onto an equal mixture {f a — 
0.5). After release and 17 ms time of flight, \b) atoms appear 
in | ±2). (B) Growth of population in |6, ±2) as determined by 
blasting away the \a) atoms after an evolution time r. Each 
data point is averaged over 6 runs. 



vant for adiabatically ramped- up, state- selective optical 
lattices. For such a system, interspecies interactions can 
be expected to give rise to diffraction effects qualitatively 
similar to those due to FWM. The density profile of the 
| a) component gets spatially modulated by the optical 
lattice, thus forming an "atomic lattice" that, in turn, 
should modulate the density of the \b) component, lead- 
ing to diffraction peaks at ±2^/cl immediately after re- 
lease. However, we find that, at least as long as both 
components are in the superfluid state, FWM is by far 
the dominant mechanism for the emergence of recoiling 
| b) atoms, caused by the projection of the \a) component 
into plane- wave momentum modes after release. We note 
that the mismatch between the dispersion relations for 
| a) and \b) atoms suppresses FWM while the lattice is 
on. For f a = 0.5 and V a = 6 Er (V& = 0), as shown 
in Fig. [4] (A), we measure a relative population of up to 
1.5 % in each of the \b, ±2) states. Assuming this to be 
caused by a density modulation would require an atomic 
lattice modulation depth of 2 E^ more than the chem- 
ical potential of the condensate in the lattice. A blast 
measurement as discussed above shows that the popula- 
tion in the observed peaks slowly grows only after release 
from the lattice [Fig.[4](B)], indicating that the peaks are 
indeed caused by FWM. 

We have performed analogous experiments for differ- 
ent final lattice depths V a , and with additional, state- 
independent lattices along the x and y directions. These 
results will be presented in detail in a future publication. 
In brief, we observe similar FWM effects along the state- 
dependent axis (cf. also [15], Fig. 8); however, we find 
that no FWM peaks are produced when the \a) atoms 
are in the Mott regime. This is consistent with the no- 
tion that FWM as described relies on the existence of a 
well-defined macroscopic phase and thus bears the poten- 
tial to be used as a sensitive probe of phase coherence. 

To summarize, we have demonstrated collinear four- 



wave mixing in a two-component mixture of bosonic 
atoms, and find excellent agreement with a simple the- 
oretical model. Our work is of relevance both in the 
context of quantum atom optics, and for experimental 
studies of bosonic mixtures in optical lattices. 
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